Remnant superfluid collective phase oscillations in 
the normal state of systems with resonant pairing 
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The signature of superfluidity in bosonic systems is a sound wave-like spectrum of the single 
particle excitations which in the case of strong interactions is roughly temperature independent. In 
fermionic systems, where fermion pairing arises as a resonance phenomenon between free fermions 
and paired fermionic states (examples are: the atomic gases of ®Li or ^''K controlled by a Feshbach 
resonance, polaronic systems in the intermediary coupling regime, d-wave hole pairing in the strongly 
correlated Hubbard system), remnants of such superfluid characteristics are expected to be visible 
in the normal state. The single particle excitations maintain there a sound wave like structure for 
wave vectors above a certain qmin{T) where they practically coincide there with the spectrum of 
the superfluid phase for T <Tc. Upon approaching the transition from above this region in q-space 
extends down to small momenta, except for a narrow region around q = Q where such modes change 
into damped free particle like excitations. 

PACS numbers: 03.75.Kk, 03.75.Ss, 74.20.Mn 



I. INTRODUCTION 

Approaching the transition to a superconducting or 
superfluid state from above, one can (under certain con- 
ditions) observe incipient macroscopic features which are 
caused by the emergence of an order parameter. In clas- 
sical superconductors such features, related to spatial 
order parameter fluctuations, are restricted to only an 
extremely narrow temperature region around the super- 
conducting critical temperature Tc, and in practice are 
hard to detect at all. Such fluctuations however are vis- 
ible in systems with real space pairing or, more gener- 
ally, when the overlap between the pair wave functions is 
small and we are in the crossover regime between a BCS 
type superfluidity of Cooperons and a superfluid phase of 
tightly bound Fermions which behave as bosons. Rem- 
nants of superfluidity, sometimes termed localized super- 
fluidity, above have been observed-'- in form of finite 
range phase correlations in purely bosonic systems such 
as liquid ^He in porous media of vicors and aerogels, with 
a characteristic disorder and confinement. Similar fea- 
tures have been seen for Fermionic systems such as '^He 
in aerogelsS and superconducting hetero structures^.. So- 
lution to the theoretical questions raised in this connec- 
tion lies in a formulation capable of describing on equal 
footing a BCS type superconductivity in a system of 
weakly coupled fermions and a Bose Einstein conden- 
sation (BEC) of strongly bound fermion pairs. Early 
attempts to do that go go back to the work of Leggett4 
and Nozieres and Schmitt-Rink^. and rely on cross-over 
scenarios where electron pairing is given by some unspec- 
ified effective attraction between them. 

Fermionic systems where the binding between 
Fermions comes about from an exchange interaction be- 
tween free itinerant Fermions and two-Fermion bound 
states, present a different scenario to examine the cross- 
over regime between a BCS type superfluidity and a con- 



densed states of tightly bound pairs. Such systems have 
moreover the advantage that sometimes, in real systems, 
the cross-over can be tuned experimentally. An exam- 
ple for such scenarios are Many Polaron systems in the 
intermediary coupling regime where free itinerant elec- 
trons engage in a resonant scattering process with weakly 
bound bipolaronic states when their respective energy 
difference is small^. This leads to long lived electron pairs 
which ultimately can condensate. An other example, now 
widely studied in the literature in connection with their 
condensation^, are gases of fermionic atomic (such as ^Li 
and *'^K atoms) which can be brought into such resonant 
fermionic pair states via a so called Feshbach resonance 
mechanism— which involves hyperfine spin-flip processes 
between the nuclear and the electronic spins of the atoms 
together with their molecular counterparts. Finally, also 
in the highly debated scenarios for the high tempera- 
ture superconductors (HTSC) resonant pairing between 
d-wave holes has been invoked. There, it has been sug- 
gested that such pairing arises from an exchange between 
itinerant holes and bound hole pairs in plaquette RVB 
states on flnite clustersiS. 

In all those systems resonant pairing leads to long lived 
electron pairs which ultimately are driven into a super- 
fluid phase. Furthermore such systems are characterized 
by a strongly interdependent dynamics of single- and two- 
particle excitations which, upon approaching and passing 
through the superconducting phase transition, simulta- 
neously undergo qualitative changes. Thus, the opening 
of a pseudogap in the single particle spectrum, when Tc 
is approached from above, occurs concomitantly with a 
change-over from single particle fermionic transport to 
one ensured by bosonic molecular entitiesi^. The ob- 
served transient Meissner effecfl^ and a Nernst effect^^ in 
the normal phase in HTSC can be considered to be signa- 
tures of that. In the atomic gases the physics is more in- 
volved because of the strong inhomogeneous character of 
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those trapped gases, leading to radial and axial breathing 
modeai^ instead of the usual sound-wave like excitation 
spectrum known in translational invariant homogeneous 
superfluids. Nevertheless, corresponding manifestation 
of superfluid fluctuations in the normal state should also 
be expected in those systems. 

In this paper we shall analyze the molecular (and/or 
fermion-pair) excitation spectrum of such a general class 
of systems which can be described in terms of resonating 
pairs of fermions and discuss how, on a finite length scale, 
superfluid phase fluctuations can emerge upon approach- 
ing Tc from above. We restrict ourselves here to the study 
of homogeneous systems, leaving the more complex struc- 
tures to be expected for remnant collective modes in in- 
homogeneous atomic gases in optical traps to a future 
work. The simplest, and generally adopted approach to 
study such systems, is on the basis of a phenomenological 
Boson-Fermion model. 

Pairing in such a model can be viewed as the An- 
dreev scattering processes between itinerant carriers and 
bosonic bound pairs on small clusters. One then is faced 
on one hand with local intra-cluster phase correlations 
between pairs of itinerant fermions and bosonic bound 
Fermion pairs and on the other hand with non-local inter- 
cluster phase correlationsi^ii^. The first ones play the 
role of local density fluctuations and the second ones 
of effective inter-site Josephson coupling. This physics, 
which is an intrinsic ingredient of the various representa- 
tive examples cited above and which are effectively real- 
ized in nature, is qualitatively different from that of the 
standard cross-over scenarios based on effective attrac- 
tive inter-particle interactions. It leads to novel features 
such as superfluid-insulator transitions, and lets one en- 
visage the possibility of normal state bose metals and 
exotic elementary as well as collective excitations which 
remain to be fully exploredi^. 

II. THE MODEL 

The following boson fermion model (BFM) Hamilto- 
nian for resonant pairing 

H ^ ^ £kc[^Ck<x + w XI (^qCq-kiCicT + h.c.) 

k,cr k,q 

+ X(£;q + 2^.)fet5q. (1) 
q 

is currently employed in studies of the above mentioned 
systems. The operators cj^^ (cko-) correspond, according 
to the physical system we are studying, to the creation 
(annihilation) of either free electrons, or free itinerant 
holons or fermionic atoms in one of two possible hyperfine 
configurations, denoted symbolically by u =t and a =[. 
The energy et of those fermions is measured with respect 
to the chemical potential /i. Correspondingly, b]^ (&q) re- 
fer to bound diatomic molecules of bosonic character (ei- 
ther localized bipolarons, or bound hole pairs on plaque- 



tte RVB states or weakly bound pairs of atoms in a triplet 
configuration), having an energy _Eq being measured with 
respect to 2/i. The parameter 2v — Eq^Q — Sskp (where 
ki? is the Fermi momentum), denotes the difference in 
energy of the weakly bound fermion pairs and the single 
fermion scattering states. If ly is small, pairing will be in- 
troduced among the uncorrelated fermions via resonance 
scattering, tantamount to a boson-fermion pair exchange 
with coupling strength v. Tuning the value of v, one can 
cover the whole regime between Cooper pairs and locally 
bound pairs and their corresponding condensed phases. 

Such a BFM has been introduced originally in solid 
state theory many years ago, in an attempt to describe 
the situation of intermediary electron - lattice couplingi^ 
and has been intensively studied over the last decade, 
mainly in connection with the pseudogap phenomenon in 
the HTSC. As shown recentl}*i2i, this model does indeed 
capture the resonant-type scattering between fermions 
due to the Feshbach mechanism and has been widely 
studied in connection with several issues of the atomic 
gas superfluidity^^. 

Our main objective here is to study the two-fermion 
dynamical correlation functions when the detuning v 
from the resonance is small, thus putting ourselves 
in the center of the cross-over regime between a 
superfluid ground state of BCS characteristics and 
one corresponding to tightly bound fermion pairs of 
bosonic character. The Green's function describing the 
fermion pairs G^""^"^ {ci, uj) is related to the single par- 
ticle boson propagator via G^(q, lj) — G^(q, lj) + 
z;2G^(q,c.)Gf^-(q,c.)G^(q,c.), where G^(q,c.) = 
[uj — Eq — 2zy] . This implies that the excitation ener- 
gies of the bound molecules and fermionic diatomic pairs 
are identical. Only the spectral weights differ as can be 
seen from the relation between their spectral functions, 
i.e., AP°*'^(q,u;) = y-^ {uj - E^ - 2iyf A^(q,w). It is 
thus sufficient to determine one of these functions in or- 
der to derive the excitation spectra for both. 



III. THE PROCEDURE 

The inter-dependence between the single and two- 
particle correlations requires to treat them on equal foot- 
ings. For that purpose we employ a continuous renor- 
malization group procedure^^ which, through a set of in- 
finitesimal canonical transformations, reduces the initial 
Hamiltonian ^ to an essentially diagonalizable form, 
containing the relevant physics which we want to de- 
scribe, plus additional terms which can be treated as 
small perturbations. Contrary to standard renormaliza- 
tion group techniques, where one integrates out the high 
energy states and subsequently derives an effective low 
energy Hamiltonian, in this method both, the high and 
low energy sectors, are renormalized and kept throughout 
the whole transformation process. 

The specific construction of such a procedure for 
the BFM was given previouslji^S,, where also the single 
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particle spectrum of the fermionic atoms, pointing to 
Bogoliubov-like excitations below as well as above T^. was 
studied^. We apply here this procedure for the study of 
the boson spectral function. In the course of diagonaliz- 
ing the Hamiltonian, the boson operators evolve toward 
a form given by 

6q = ^qfeq+-^^6q,kCkxCq_kT, ^q=(^q)^- (2) 

The two complex coefficients appearing in are cal- 
culated in the limit of the convergence of the renor- 
malization flow procedure limj^oo -4q(/) = Aq and 

linij^oo 'Sq^klO = '^q.kj where I denotes the continuous 
flow parameter. We base ourselves on the general rela- 
tions which describe the evolution of operators^ 

dO{l)/dl = [7il),0{l)] (3) 

with r] being judiciously choseniS as 



k,p 



and ak,p(0 = [ek(0 + ep(0 - ^^k+pCO] VKpi^)^- CoefH- 
cients Aq{l) and -Bq_k(0 satisfy then the renormalization 
equations 



dAqjl) 
dl 

dBq,u{l) 



^ -^E"k.q-k(0/k,q-kSq,k(0, (5) 
k 



dl 



ak,q-k(0-^q(0: 



(6) 



with the initial conditions Aq{Q) = 1, i3q.k(0) = and 

fKP = - Ki - ripr 

This procedure leads finally to the following form of 

the spectral function for the bosonic molecules 



(7) 



A^(q,c^) = \Aq\'S{LU~E, 



/k.q-kl^q.kPl^ (W - £k - £q-k) • 



The first term of |7J describes long-lived quasi-particles 
with the renormalized energy and whose spectral 

weight is |-4qp. The second term describes the inco- 
herent background extending over the region determined 
by the renormalized fermion energies Syr^. From equa- 
tions (|5I6|) we derive the following sum rule |-4qP -I- 
;^ Sk l'^q^kP/k,q-k = 1 which correctly preserves the 
total spectral weight f^^diLjA^{q,uj) = ([foq,foq]) = 1- 



IV. THE PAIR EXCITATION SPECTRUM 
BELOW Tc 

At a certain critical temperature the static pair sus- 
ceptibility Xlk.p lo '^^e™(4^(T)c|j_^^(r)cq_pxCpt)|,^^o 




the long-lived boson and/or fermion pair excitation spectrum 
in the superfluid state at T = 0. The shaded regions show 
the incoherent background, which is energetically separated 
from the collective excitation branch. At q = the incoher- 
ent background exists for energies larger than 2vxb (twice 
the value of the single particle fermion gap). We used the fol- 
lowing dispersions Ck = — ■jCos(a|k|) and Eq = — -jC0s(a|q|) 
such that the mass ms = 2mF and the potential v = O.ID. 
We further set the lattice constant a = 1 and use the band- 
width D as a unit for energies. 



becomes divergent for q = and, due to the Thou- 
less criterion, the system undergoes a phase transition 
to a superfluid state. For T < there appear two 
order parameters which are proportional to each other: 
Xf = (c-kj,Ck|) for the fermions and xb = (^q=o) for the 
bosons (atom molecules). 

Near the Fermi energy, the single parti- 
cle fermionic excitations become gaped: Ek — 
sgn {ek} V (^k)^ + {vXbY- In consequence, no fermionic 
states, neither coherent nor incoherent, exist within 
the energy window < vx^^- This simultaneously 
affects the incoherent part of bosonic spectrum, as can 
be seen from (UJ. For the long wavelength limit q ^ 
the incoherent background is pushed up to energies 
\uj\ > 2vxb and thus permits long-lived excitations, 
which correspond to collective modes, known as first 
sound for interacting bosonic systems in the superfluid 
state (see Fig.l). The temperature dependence of these 
modes has previously been studied for this BFMM in 
the superfluid phase within a framework of the dielectric 
formalism with use of the Ward identities, currently 
employed in the theory of interacting bose gases. A 
behavior similar to that of the strong coupling limit of 
interacting bose gases^^ was found, showing a sound 
velocity being little dependent on temperature as one 
traverses the superfluid transition, but whose spectral 
weight in the boson single particle spectral function 
disappears upon approaching Tc- 

Such sound wave-like modes are not realized in charged 
superconducting systems because of the long range 
Coulomb interaction which pushes them up to the gen- 
erally huge plasma frequencjjS^. For electrically neutral 
atoms, such as the trapped atomic gases, this is no longer 
the case and hence one can realistically expect collective 
sound wave-like modes, although appropriately modified 
due to the inhomogeneous structure of the gas densitjii^. 
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FIG. 2: Comparison of the dispersion Eq of the coherent part 
of the boson spectral function at temperatures correspond- 
ing to the superfluid (T = 0), pseudogap (0.007, 0.01) and 
the normal phase above T* (0.02). Upper panel shows the 
derivative dEq/dq of these curves. The insets contain corre- 
spondingly: the low momentum q limit of the dispersion Eq 
and the temperature dependence of chemical potential (the 
marked points correspond to four temperatures T = 0.02, 
0.01, 0.007 and chosen in this work). 



V. THE PAIR EXCITATION SPECTRUM 
ABOVE Tc 

Decreasing the temperature in the normal state below 
a certain T*(> Tc) one expects precursor pairing effects 
which show up in the single particle fermionic excitations 
spectrum in form of a pseudogap which opens up near 
the chemical potential^SiSiiS^. Above T* the low energy 
part of the pair excitations has the usual parabolic dis- 
persion. However, upon decreasing the temperature and 
approaching Tc, phase coherence gradually sets in on a fi- 
nite length and time scale, which becomes visible in form 
of a linear in q dispersion of the single boson (respectively 
Fermion pair) excitation for small q vectors, in an inter- 
val [qmin{T),qmax{T)] (scc Fig. 2). There, the deriva- 
tive of the effective Bose single particle energy spectrum 
dEq/dq shows a flat portion, which, when extrapolated 
to q = 0, practically coincides with the corresponding 
quantity in the superfluid phase at T = 0. We observe 
that, as the temperature is decreased, qmin{T) decreases 
toward zero, but always leaving a small interval in q space 
[0, qmin{T)\ where one clearly observes a free particle like 
spectrum with an effective mass which decreases as T de- 
creases. This is in accordance with an earlier study on 
this subject using selfconsistent perturbation theory?^. 
For T >T* the coherent boson mode overlaps with an in- 
coherent background in the single particle boson spectral 
function (see Fig. 3, bottom panel). Yet, upon decreas- 
ing the temperature to below T*, we observe that this 
incoherent background moves away from the position of 
the coherent contribution (upper panel of Fig. 3) which 
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FIG. 3: The boson spectral function A^{(\,ui) for the low 
energy pair excitations. The upper panel corresponds to T = 
0.007 being close to Tc in the pseudogap region T* > T > 
Tc. The bottom panel refers to the normal state T — 0.02 
(above T*). In the pseudogap phase a propagating coherent 
contribution given by the 5-function peak and an incoherent 
background, given by the shaded regions, get separated above 
some relatively small critical momentum qmin- This is no 
longer the case for T > T* . 



ensures that a linear in q branch of the boson spectrum 
is well defined in the corresponding interval of q vectors. 
This strongly suggests that remnants of the first sound 
still can exist as part of the single particle boson spec- 
trum above the superfluid phase transition for a limited 
region of wave vectors due to a persistence of superfluid 
phase correlations above Tc on a flnite length and time 
scale. 



VI. CONCLUSIONS 

We studied the qualitative changes of the excitation 
spectrum for the resonant fermion pairs which occur 
upon varying the temperature. We found that quantum 
fluctuations play a crucial role when detuning v from 
the Feshbach resonance is small. Fluctuations manifest 
themselves in the pseudogap regime T* > T > Tc- 

Far above Tc the off-diagonal long range order is not 
established. The pair excitation spectrum for small q 
vectors is then characterized by a parabolic branch (see 
Fig. 2) and overlaps with the incoherent background (see 



5 



the bottom panel of Fig. 3) such as to effectively destroy 
any bosonic quasi-particle features. 

This situation changes dramatically when the temper- 
ature drops below T* where resonant pairing sets in. 
Phase correlations start to build up on a finite spatial 
and temporal scale as the temperature decreases and ap- 
proaches Tc. The single particle fermion spectrum re- 
veals then a partial suppression of states (pseudogap) 
around the Fermi energji2Si2L2^, which is accompanied 
by qualitative changes in the pair excitation spectrum. 
Quantum fluctuations lead to emergence of the collec- 
tive sound-wave mode which above exists in a finite 
momentum interval [qminiT), qraax{T)]. Upon decreasing 
the temperature the long-lived branch of the pair spec- 
trum gradually splits off from the incoherent background 
(upper panel in Fig. 3) and spreads over a wider and 
wider momentum region, with QminiT) steadily decreas- 
ing as we approach Tc- We note however that invariably 
the linear in q dispersion changes into a damped free par- 
ticle like behavior in the close vicinity of g = 0. 

The sound-wave mode has been so far measured above 
Tc in the liquid helium by ultrasonic techniquesSl as 



well as by neutron scattering measurementai. In the 
case of the trapped fermionic atoms the corresponding 
mode is expected to be compressional density waves and, 
similar to the present study, one should expect rem- 
nants of those modes in the normal state. In principle, 
such modes can be experimentally checked by the Bragg 
spectroscopy^^. Indirect methods for detecting the col- 
lective modes which rely on measuring the magnetic sus- 
ceptibility and density-density correlation functions have 
been discussed (although only for T < Tc) in Refi^. In 
some future work we shall discuss how collective modes 
can possibly be observed in measurements of the mag- 
netic susceptibility in the pseudogap regime above Tc- 
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